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Abstract 

In this paper we survey the primary research, both theoretical and applied, in the area of Robust 
Optimization (RO). Our focus is on the computational attractiveness of RO approaches, as well as 
the modeling power and broad applicability of the methodology. In addition to surveying prominent 
theoretical results of RO, we also present some recent results linking RO to adaptable models for 
multi-stage decision-making problems. Finally, we highlight applications of RO across a wide spec- 
trum of domains, including finance, statistics, learning, and various areas of engineering. 
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1 Introduction 



Optimization affected by parameter uncertainty has long been a focus of the mathematical programming 
community. Solutions to optimization problems can exhibit remarkable sensitivity to perturbations in 
the parameters of the problem (demonstrated in compelling fashion in [15]) thus often rendering a 
computed solution highly infeasible, suboptimal, or both (in short, potentially worthless). 

In science and engineering, this is hardly a new notion. In the context of optimization, the most 
closely related field is that of Robust Control (we refer to the textbooks [T36J and |67j , and the references 
therein). While there are many high-level similarities, and indeed much of the motivation for the de- 
velopment of Robust Optimization came from the Robust Control community, Robust Optimization is 
a distinct field, focusing on traditional optimization-theoretic concepts, particularly algorithms, geom- 
etry, and tractability, in addition to modeling power and structural results which are more generically 
prevalent in the setting of robustness. 

Stochastic Optimization starts by assuming the uncertainty has a probabilistic description. This 
approach has a long and active history dating at least as far back as Dantzig's original paper [61J. We 
refer the interested reader to several textbooks ([90 1 l39 l 11211 [92] ) and the many references therein for a 
more comprehensive picture of Stochastic Optimization. 

This paper considers Robust Optimization (RO), a more recent approach to optimization under 
uncertainty, in which the uncertainty model is not stochastic, but rather deterministic and set-based. 
Instead of seeking to immunize the solution in some probabilistic sense to stochastic uncertainty, here 
the decision-maker constructs a solution that is feasible for any realization of the uncertainty in a 
given set. The motivation for this approach is twofold. First, the model of set-based uncertainty is 
interesting in its own right, and in many applications is an appropriate notion of parameter uncertainty. 
Second, computational tractability is also a primary motivation and goal. It is this latter objective that 
has largely influenced the theoretical trajectory of Robust Optimization, and, more recently, has been 
responsible for its burgeoning success in a broad variety of application areas. The work of Ben-Tal and 
Nemirovski (e.g., [131 El US]) and El Ghaoui et al. [TTJ |HD] in the late 1990s, coupled with advances 
in computing technology and the development of fast, interior point methods for convex optimization, 
particularly for semidefinite optimization (e.g., Boyd and Vandenberghe, |42j) sparked a massive flurry 
of interest in the field of Robust Optimization. 

Central issues we seek to address in this paper include tractability of robust optimization models; 
conservativeness of the RO formulation, and flexibility to apply the framework to different settings and 
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applications. We give a summary of the main issues raised, and results presented. 

1. Tractability: In general, the robust version of a tractabk£] optimization problem may not itself be 
tractable. We outline tractability results, which depend on the structure of the nominal problem as 
well as the class of uncertainty set. Many well-known classes of optimization problems, including 
LP, QCQP, SOCP, SDP, and some discrete problems as well, have a RO formulation that is 
tractable. Some care must be taken in the choice of the uncertainty set to ensure that tractability 
is preserved. 

2. Conservativeness and probability guarantees: RO constructs solutions that are deterministically 
immune to realizations of the uncertain parameters in certain sets. This approach may be the 
only reasonable alternative when the parameter uncertainty is not stochastic, or if distributional 
information is not readily available. But even if there is an underlying distribution, the tractabil- 
ity benefits of the Robust Optimization approach may make it more attractive than alternative 
approaches from Stochastic Optimization. In this case, we might ask for probabilistic guarantees 
for the robust solution that can be computed a priori, as a function of the structure and size 
of the uncertainty set. In the sequel, we show that there are several convenient, efficient, and 
well-motivated parameterizations of different classes of uncertainty sets, that provide a notion of 
a budget of uncertainty. This allows the designer a level of flexibility in choosing the tradeoff 
between robustness and performance, and also allows the ability to choose the corresponding level 
of probabilistic protection. In particular, a perhaps surprising implication is that while the robust 
optimization formulation is inherently max-min (i.e., worst-case), the solutions it produces need 
not be overly conservative, and in many cases are very similar to those produced by stochastic 
methods. 

3. Flexibility: In Section[2j we discuss a wide array of optimization classes, and also uncertainty sets, 

and consider the properties of the robust versions. In the final section of this paper, we illustrate 

the broad modeling power of Robust Optimization by presenting a wide variety of applications. 

We also give pointers to some surprising uses of robust optimization, particularly in statistics, 

1 Throughout this paper, we use the term "tractable" repeatedly. We use this as shorthand to refer to problems that 
can be reformulated into equivalent problems for which there are known solution algorithms with worst-case running 
time polynomial in a properly defined input size (see, e.g., Section 6.6 of Ben-Tal and Nemirovski 19 ). Similarly, by 
"intractable" we mean the existence of such an algorithm for general instances of the problem would imply P=NP. 
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where robust optimization is used as a tool to imbue the solution with desirable properties, like 
sparsity, stability or statistical consistency. 

The overall aim of this paper is to outline the development and main aspects of Robust Optimization, 
with an emphasis on its flexibility and structure. While the paper is organized around some of the 
main themes of robust optimization research, we attempt throughout to compare with other methods, 
particularly stochastic optimization, thus providing guidance and some intuition on when the robust 
optimization avenue may be most appropriate, and ultimately successful. 

We also refer the interested reader to the recent book of Ben-Tal, El Ghaoui and Nemirovski |19j . 
which is an excellent reference on Robust Optimization that provides more detail on specific formulation 
and tractability issues. Our goal here is to provide a more condensed, higher level summary of key 
methodological results as well as a broad array of applications that use Robust Optimization. 

A First Example 

To motivate RO and some of the modeling issues at hand, we begin with an example from portfolio 
selection. The example is a fairly standard one. We consider an investor who is attempting to allocate 
one unit of wealth among n risky assets with random return r and a risk-free asset (cash) with known 
return r/. The investor may not short-sell risky assets or borrow. His goal is to optimally trade off 
between expected return and the probability that his portfolio loses money. 

If the returns are stochastic with known distribution, the tradeoff between expected return and loss 
probability is a stochastic program. However, calculating a point on the pareto frontier is in general 
NP-hard even when the distribution of returns is discrete (Benati and Rizzi |20j). 

We will consider two different cases: one where the distribution of asset price fluctuation matches 
the empirical distribution of given historical data and hence is known exactly, and then the case where 
it only approximately matches historical data. The latter case is of considerable practical importance, 
as the distribution of new returns (after an allocation decision) often deviate significantly from past 
samples. We compare the stochastic solution to several easily solved RO-based approximations in both 
of these cases. 

The intractability of the stochastic problem arises because of the probability constraint on the loss: 



where x is the vector of allocations into the n risky assets (the decision variables). The robust opti- 
mization formulations replace this probabilistic constraint with a deterministic constraint, requiring the 
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return to be nonnegative for any realization of the returns in some given set, called the uncertainty set: 

r'x + r/(l - l'aj) > 1 Vf € 11. (1.2) 



While not explicitly specified in the robust constraint (1.2), the resulting solution has some pi oss . As 
a rough rule, the bigger the set 1Z, the lower the objective function (since there are more constraints 
to satisfy) , and the smaller the loss probability pi OS s • Central themes in robust optimization are under- 
standing how to structure the uncertainty set 1Z so that the resulting problem is tractable and favorably 
trades off expected return with loss probability pi oss ■ Section [2] is devoted to the tractability of different 
types of uncertainty sets. Section [3] focuses on obtaining a priori probabilistic guarantees given differ- 
ent uncertainty sets. Here, we consider three types of uncertainty sets, all defined with a parameter to 
control "size" so that we can explore the resulting tradeoff of return, and pi oss : 

TZ Q (j) = {r : (r-r)'^- l {f-r) < j 2 } , 

Tl D {T) = Ir :3u£Rl s.t. r i =f i + (r i -f i )u i , Ui<l,^2ui<r\, 



i=l 

N 



K T (a) = if : 3g€R? s.t.f = ^ gi rS l'q = l, gi < N ^_ a y » = 1, ■ ■ • j ■ 

The set TZ Q (j) is a quadratic or ellipsoidal uncertainty set: this set considers all returns within a 
radius of 7 from the mean return vector, where the ellipsoid is tilted by the covariance. When 7 = 0, this 
set is just the singleton {r }. The set 1Z D (T) (D for "D-norm" model considered in Section [2]) considers 
all returns such that each component of the return is in the interval fi], with the restriction that 
the total weight of deviation from fj, summed across all assets, may be no more than T. When T = 0, 
this set is the singleton {f}; at the other extreme, when V = n, returns in the range [r^i] for all 
assets are considered. Finally, 7Z T (k) is the "tail" uncertainty set, and considers the convex hull of all 
possible N(l — a) point averages of the N returns. When a = 0, this set is the singleton {f }. When 
a = (N — 1)/N, this set is the convex hull of all iV returns. 

To illustrate the use of these formulations, consider n = 10 risky assets based on iV = 300 past 
market returns. The assets are a collection of equity and debt indices, and the return observations are 
monthly from a data set starting in 1981. For each of the three uncertainty RO formulations, we solve 
200 problems, each maximizing expected return subject to feasibility and the robust constraint at one 
of 200 different values of their defining parameter 7, T, or a. In total, we solve 600 RO formulations. 
For comparison, we also formulate the problem of minimizing probability of loss subject to an expected 
return constraint as a stochastic program (which can be formulated as a mixed integer program), and 
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Figure 1: (L): Expected return- loss probability frontier for RO-based formulations and exact stochastic 
formulation; numbers are time (sec.) for solving each stochastic program. (C/R): Frontier for model 
with random perturbations bounded by 1% (C) and 2% (R). 

solve 8 versions of this problem, each corresponding to one of 8 different expected return levels. The 
computations are performed using the MOSEK optimization toolbox in Matlab on a laptop computer 
with a 2.13GHZ processor and 2GB of RAM. 

The results are shown in Figure [T] On the left, we see the frontier for the three RO-based formulations 
as well as the performance of the exact formulation (at the 8 return levels). The numbers indicate the 
time in seconds to solve the stochastic program in each case. 

The stochastic model is designed for the nominal case, so we expect it to outperform the three 
RO-based formulations. However, even under this model, the gap from the 'RP and 1Z T RO frontiers is 
small: in several of the cases, the difference in performance is almost negligible. The largest improvement 
offered by the stochastic formulation is around a 2 — 3% decrease in loss probability. Here, the solutions 
from the TZr model do not fare as well; though there is a range in which its performance is comparable to 
the other two RO-based models, typically its allocations appear to be conservative. In general, solving 
the stochastic formulation exactly is difficult, which is not surprising given its NP-hardness. Though a 
few of the instances at extreme return levels are solved in only a few seconds, several of the instances 
require well over an hour to solve, and the worst case requires over 2.1 hours to solve. The total time 
to solve these 8 instances is about 5.2 hours; by contrast, solving the 600 RO-based instances takes a 
bit under 10 minutes in total, or about one second per instance. 

On the center and right parts of Figure [T] are results for the computed portfolios under the same 
return model but with random perturbations. Specifically, we perturb each of the N x n return values 
by a random number uniformly distributed on [.99, 1.01] in the middle figure and [.98, 1.02] in the right 
figure. At the 1% perturbation level, the gap in performance between the models is reduced, and there 
are regions in which each of the models are best as well as worst. The model based on 1Z D is least 
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affected by the perturbation; its frontier is essentially unchanged. The models based on and TZ T 
are more significantly affected, perhaps with the effect on 1Z T being a bit more pronounced. Finally, 
the stochastic formulation's solutions are the most sensitive of the bunch: though the SP solution is a 
winner in one of the 8 cases, it is worse off than the others in several of the other cases, and the increase 
in loss probability from the original model is as large as 5 — 6% for the SP solutions. 

At the 2% level, the results are even more pronounced: here, the SP solutions are always outper- 
formed by one of the robust approaches, and the solutions based on 1Z D are relatively unaffected by 
the noise. The other two robust approaches are substantially affected, but nonetheless still win out in 
some parts of the frontier. When noise is introduced, it does not appear that the exact solutions confer 
much of an advantage and, in fact, may perform considerably worse. Though this is only one random 
trial, such results are typical. 

There are several points of discussion here. First is the issue of complexity. The RO-based models 
are all fairly easy to solve here, though they themselves have complexities that scale differently. The 1Z® 
model may be formulated as a second-order cone program (SOCP); both the TZ D and the TZ T models 
may be formulated as an LP. Meanwhile, the exact stochastic model is an NP-hard mixed integer 
program. Under the original model, it is clearly much easier to solve these RO-based models than the 
exact formulation. In a problem with financial data, it is easy to imagine having thousands of return 
samples. Whereas the RO formulations can still be solved quickly in such cases, solving the exact SP 
could be hopeless. 

A second issue is the ability of solution methods to cope with deviations in the underlying model 
(or "model uncertainty"). The RO-based formulations themselves are different in this regard. Here, 
the 1Z D approach focuses on the worst-case returns on a subset of the assets, the 1Z® approach focuses 
on the first two moments of the returns, and the TZ T approach focuses on averages over the lower tail 
of the distribution. Though all of these are somehow "robust," 1Z D is the "most robust" of the three; 
indeed, we also implemented perturbations at the 5% level and found its frontier is relatively unchanged, 
while the other three frontiers are severely distorted. Intuitively, we would expect models that are more 
robust will fare better in situations with new or altered data; indeed, we will later touch upon some 
work that shows that there are intimate connections between the robustness of a model and its ability 
to generalize in a statistical learning sense. 

This idea - that Robust Optimization is useful in dealing with erroneous or noise-corrupted data 
- seems relatively well understood by the optimization community (those who build, study, and solve 
optimization models) at-large. In fact, we would guess that many figure this to be the raison d'etre 
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for Robust Optimization. The final point that we would like to make is that, while dealing with 
perturbations is one virtue of the approach, RO is also more broadly of use as a computationally viable 
way to handle uncertainty in models that are on their own quite difficult to solve, as illustrated here. 

In this example, even if we are absolutely set on the original model, it is hard to solve exactly. 
Nonetheless, two of the RO-based approaches perform well and are not far from optimal under the 
nominal model. In addition, they may be computed orders of magnitude faster than the exact solution. 
Of course, we also see that the user needs to have some understanding of the structure of the uncertainty 
set in order to intelligently use RO techniques: the approach with 1Z D , though somewhat conservative 
in the original model, is quite resistant to perturbations of the model. 

In short, RO provides a set of tools that may be useful in dealing with different types of uncertainties 
- both the "model error" or "noisy data" type as well as complex, stochastic descriptions of uncertainty 
in an explicit model - in a computationally manageable way. Like any approach, however, there are 
tradeoffs, both in terms of performance issues and in terms of problem complexity. Understanding and 
managing these tradeoffs requires expertise. The goal of this paper, first and foremost, is to describe 
some of this landscape for RO. This includes detailing what types of RO formulations may be efficiently 
solved in large scale, as well what connections various RO formulations have to perhaps more widely 
known methods. The second goal of this paper is to then highlight an array of application domains for 
which some of these techniques have been useful. 

2 Structure and tractability results 

In this section, we outline several of the structural properties, and detail some tractability results of 
Robust Optimization. We also show how the notion of a budget of uncertainty enters into several 
different uncertainty set formulations. 

2.1 Robust Optimization 

Given an objective fo(x) to optimize, subject to constraints fi(x,Ui) < with uncertain parameters, 
{ui}, the general Robust Optimization formulation is: 

minimize fo( x ) 

subject to fi{x, < 0, V u, 6 Ui, i = 1, . . . ,m. (2.3) 
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Here x G M n is a vector of decision variables, fo, fi ■ K n — > K are functions, and the uncertainty 
parameters £ M fc are assumed to take arbitrary values in the uncertainty sets Ui C which, for 



our purposes, will always be closed. The goal of (2.3) is to compute minimum cost solutions x* among 



all those solutions which are feasible for all realizations of the disturbances u» within Ui. Thus, if some 



of the lAi are continuous sets, (2.3), as stated, has an infinite number of constraints. Intuitively, this 
problem offers some measure of feasibility protection for optimization problems containing parameters 
which are not known exactly. 



It is worthwhile to notice the following, straightforward facts about the problem statement of (2.3): 



The fact that the objective function is unaffected by parameter uncertainty is without loss of 
generality; we may always introduce an auxiliary variable, call it t, and minimize t subject to the 
additional constraint max fo(x,uo) < t. 

It is also without loss of generality to assume that the uncertainty set U has the form U = 
U\ x . . . x lira- If we have a single uncertainty set U for which we require (ui, . . . , u m ) £ U, then 



the constraint- wise feasibility requirement implies an equivalent problem is (2.3) with the U L taken 
as the projection of U along the corresponding dimensions (see Ben-Tal and Nemirovski, |14j for 
more on this). 

Constraints without uncertainty are also captured in this framework by assuming the correspond- 
ing Ui to be singletons. 



• Problem (2.3) also contains the instances when the decision or disturbance vectors are contained 
in more general vector spaces than W 1 or R k (e.g., 8 n in the case of semidefinite optimization) 
with the definitions modified accordingly. 

Robust Optimization is distinctly different than sensitivity analysis, which is typically applied as a 
post-optimization tool for quantifying the change in cost of the associated optimal solution with small 
perturbations in the underlying problem data. Here, our goal is to compute fixed solutions that ensure 
feasibility independent of the data. In other words, such solutions have a priori ensured feasibility when 
the problem parameters vary within the prescribed uncertainty set, which may be large. We refer the 
reader to some of the standard optimization literature (e.g., Bertsimas and Tsitsiklis, [37], Boyd and 
Vandenberghe, [33]) and works on perturbation theory (e.g., Freund, [75], Renegar, [123J) for more on 
sensitivity analysis. 
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It is not at all clear when (2.3) is efficiently solvable. One might imagine that the addition of robust- 
ness to a general optimization problem comes at the expense of significantly increased computational 
complexity. It turns out that this is true: the robust counterpart to an arbitrary convex optimization 
problem is in general intractable (|13j: some approximation results for robust convex problems with 
a conic structure are discussed in [33]). Despite this, there are many robust problems that may be 
handled in a tractable manner, and much of the literature has focused on specifying classes of functions 
fi, coupled with the types of uncertainty sets Ui, that yield tractable robust counterparts. If we define 
the robust feasible set to be 

X(U) = {x | fi(x, Ui ) < V Ui eU h i = l,...,m}, (2.4) 

then for the most partj^] tractability is tantamount to X(U) being convex in x, with an efficiently 
computable separation test. More precisely, in the next section we show that this is related to the 
structure of a particular subproblem. We now present an abridged taxonomy of some of the main 
results related to this issue. 

2.2 Robust linear optimization 

The robust counterpart of a linear optimization problem is written, without loss of generality, as 

minimize c T x 

subject to Ax < b, V a\ E U\, . . . , a m E U m , (2-5) 

where aj represents the i th row of the uncertain matrix A, and takes values in the uncertainty set 
Ui C W 1 . Then, ajx < bi, Va.; E Ui, if and only if max| a . g ^.} aj x < bi, Vi. We refer to this as the 
subproblem which must be solved. Ben-Tal and Nemirovski [14] show that the robust LP is essentially 
always tractable for most practical uncertainty sets of interest. Of course, the resulting robust problem 
may no longer be an LP. We now provide some more detailed examples. 

Ellipsoidal Uncertainty: Ben-Tal and Nemirovski [H], as well as El Ghaoui et al. [77J [80] . 
consider ellipsoidal uncertainty sets. Controlling the size of these ellipsoidal sets, as in the theorem 
below, has the interpretation of a budget of uncertainty that the decision-maker selects in order to 
easily trade off robustness and performance. 



2 i.e., subject to a Slater condition. 
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Theorem 1. (Ben-Tal and Nemirovski, WM) LetU be "ellipsoidal," i.e., 

U = U{U, Q) = {U(u) | ||Qu|| < p} , 



where u — > Tl(u) is an affine embedding of K into M mxn and Q G K . T/ien Problem (2.5) is 



equivalent to a second-order cone program (SOCP). Explicitly, if we have the uncertain optimization 

minimize cx 

subject to aiX < b{, Vaj £ Ui, Vi = 1, . . . , m, 
where the uncertainty set is given as: 

U = {(ai,...,a m ) : ai = a° + Ai-Ui, i = l,...,m, \\u\\ 2 < p}, 

((K denotes the nominal value) then the robust counterpart is: 

mininize cJ x 

subject to a®x < b{ — p\\A{x\\2, Vi = 1, . . . ,m. 

The intuition is as follows: for the case of ellipsoidal uncertainty, the subproblem max| a . g ^.| ajx < 
bi, Vi, is an optimization over a quadratic constraint. The dual, therefore, involves quadratic functions, 
which leads to the resulting SOCP. 



Polyhedral Uncertainty: Polyhedral uncertainty can be viewed as a special case of ellipsoidal un- 
certainty |14j . When U is polyhedral, the subproblem becomes linear, and the robust counterpart is 
equivalent to a linear optimization problem. To illustrate this, consider the problem: 

min : c T x 

s.t. : max {D . a .< di} aj x <bi, i = l,...,m. 
The dual of the subproblem (recall that x is not a variable of optimization in the inner max) becomes: 

min : pj di 



max : ajx 
s.t 



Didi < d, 



< — > 



s.t. 



Pi > o. 



x 



and therefore the robust linear optimization now becomes: 

min : c T x 

s.t. : pjdi <bi, i = 



Pi>0, 



x. 



, 771 

, m 



1, . . . ,m. 
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Thus the size of such problems grows polynomially in the size of the nominal problem and the dimen- 
sions of the uncertainty set. 

Cardinality Constrained Uncertainty: Bertsimas and Sim ( |34j ) define a family of polyhedral 
uncertainty sets that encode a budget of uncertainty in terms of cardinality constraints: the number 
of parameters of the problem that are allowed to vary from their nominal values. The uncertainty set 
TZP from our introductory example, is an instance of this. More generally, given an uncertain matrix, 
A = (aij), suppose each component lies in [aij — aij,a%j + a^]. Rather than protect against the 
case when every parameter can deviate, as in the original model of Soyster ([128]), we allow at most 
Ti coefficients of row i to deviate. Thus the positive number Tj denotes the budget of uncertainty for 
the i th constraint, and just as with the ellipsoidal sizing, controls the trade-off between the optimality 
of the solution, and its robustness to parameter perturbation^ Given values ri,...,r m , the robust 
formulation becomes: 

min : cx 

s.t. : J2j atjXj + maxi^cj, : |Sj|=r 4 } E je s t *aVi ^ 6 * l<i<m 

-Uj < Xj <Uj 1 < j < n ( 2 -6) 

I < x < u 

y>o. 

Because of the set-selection in the inner maximization, this problem is nonconvex. However, one can 
show that the natural convex relaxation is exact. Thus, relaxing and taking the dual of the inner 
maximization problem, one can show that the above is equivalent to the following linear formulation, 
and therefore is tractable (and moreover is a linear optimization problem): 

max : c T x 

s.t. : Yuj a ij x j + z i^i + J2jPij ^ bi Vi 
H +Pij > aijUj Vi, j 

H i < x 3 < Vi v i 

I < X < u 

V > o 
y > o. 

Norm Uncertainty: Bertsimas et al. [31] show that robust linear optimization problems with uncer- 
3 For the full details see [34]. 
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tainty sets described by more general norms lead to convex problems with constraints related to the 
dual norm. Here we use the notation vec(A) to denote the vector formed by concatenating all of the 
rows of the matrix A. 

Theorem 2. (Bertsimas et al., fSlf) With the uncertainty set 

U = {A | ||M(vec(A) - vec(A))|| < A}, 



where M is an invertible matrix, A is any constant matrix, and \\ ■ \\ is any norm, Problem (2.5) is 
equivalent to the problem 



minimize c T x 



subject to A { x + A||(M' ) a^H* < b u i = l,...,m, 

where X{ G R( m ' n ) xl is a vector that contains x G W 1 in entries (i — 1) • n + 1 through i ■ n and 
everywhere else, and \\ ■ \\* is the corresponding dual norm of \\ ■ \\. 

Thus the norm-based model shown in Theorem [2] yields an equivalent problem with corresponding 
dual norm constraints. In particular, the l\ and l^ norms result in linear optimization problems, and 
the h norm results in a second-order cone problem. 

In short, for many choices of the uncertainty set, robust linear optimization problems are tractable. 

2.3 Robust quadratic optimization 

Quadratically constrained quadratic programs (QCQP) have defining functions fi{x,u{) of the form 

fi{x,Ui) = \\Aix\\ 2 + bjx + a. 

Second order cone programs (SOCPs) have 

fi(x,Ui) = \\Aix + bj|| - cjx - di. 

For both classes, if the uncertainty set IA is a single ellipsoid (called simple ellipsoidal uncertainty) the 
robust counterpart is a semidefinite optimization problem (SDP). If IA is polyhedral or the intersection 
of ellipsoids, the robust counterpart is NP-hard (Ben-Tal and Nemirovski, |13 } 114 4 135]). 

Following [18J, we illustrate here only how to obtain the explicit reformulation of a robust quadratic 
constraint, subject to simple ellipsoidal uncertainty. Consider the quadratic constraint 

x T A T Ax<2b T x + c, V(A,b,c)eU, (2.7) 
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where the uncertainty set IA is an ellipsoid about a nominal point (A°,b°,c°): 



U 



A 



UA,b,c) := (A°,b ,c ) + J2 Ul (A l ,b l ,c l ) : \\u\\ 2 < lj . 



< 0. 



As in the previous section, a vector x is feasible for the robust constraint (2.7) if and only if it is feasible 
for the constraint: 

max : x T A T Ax — 2b T x — c 
s.t. : (A,b,c) G U 

This is the maximization of a convex quadratic objective (when the variable is the matrix A, x T A T Ax 
is quadratic and convex in A since xx T is always semidefinite) subject to a single quadratic constraint. 
It is well-known that while this problem is not convex (we are maximizing a convex quadratic) it 
nonetheless enjoys a hidden convexity property (for an early reference, see Brickman |44j ) that allows 
it to be reformulated as a (convex) semidefinite optimization problem. This is related to the so-called 
5-lemma (or S'-procedure) in control (e.g., Boyd et al. [3D], Polik and Terlaky |119j ): 

The S'-lemma essentially gives the boundary between what we can solve exactly, and where solving 
the subproblem becomes difficult. If the uncertainty set is an intersection of ellipsoids or polyhedral, 
then exact solution of the subproblem is NP-hard|^] 

Taking the dual of the SDP resulting from the S'-lemma, we have an exact, convex reformulation of 
the subproblem in the RO problem. 



Theorem 3. Given a vector x, it is feasible to the robust constraint (2.1) if and only if there exists a 
scalar r£l such that the following matrix inequality holds: 

( „0 , o™T^0 _ 1„1 i „Tj.l „L i „TuL /^O^T \ 



c° + 2x T b°-T 


lc 1 + x T b 1 ■ 


■ c L + x T b L 


(A°x) T 


lc 1 +x T b 1 


T 




(A 1 x) T 


±c L + x T b L 




T 


(A L x) T 


A°x 


A 1 x 


A L x 


I 



y o. 



2.4 Robust Semidefinite Optimization 

With ellipsoidal uncertainty sets, robust counterparts of semidefinite optimization problems are, in 
general, NP-hard (Ben-Tal and Nemirovski, [13] , Ben-Tal et al. [H]). Similar negative results hold even 
in the case of polyhedral uncertainty sets (Nemirovski, |110j ). One exception (Boyd et al. [40j) is when 



Nevertheless, there are some approximation results available: [18] 
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the uncertainty set is represented as unstructured norm-bounded uncertainty. Such uncertainty takes 
the form 

A (x) + L'{x)CR(x) + R(x)CL'(x), 

where £ is a matrix with norm satisfying HCII22 < 1> L and R are affine in the decision variables x, and 
at least one of L or R is independent of x. 

In the general case, however, robust SDP is an intractable problem. Computing approximate solu- 
tions, i.e., solutions that are robust feasible but not robust optimal to robust semidefinite optimization 
problems has, as a consequence, received considerable attention (e.g., [80], [HUTS], and [35] ). These 
methods provide bounds by developing inner approximations of the feasible set. The goodness of the 
approximation is based on a measure of how close the inner approximation is to the true feasible set. 
Precisely, the measure for this is: 

p(AR : R) = inf {p > 1 | X(AR) D X{U{p))} , 

where -Y(AR) is the feasible set of the approximate robust problem and X(U(p)) is the feasible set of 
the original robust SDP with the uncertainty set "inflated" by a factor of p. When the uncertainty 
set has "structured norm bounded" form, Ben-Tal and Nemirovski |17] develop an inner approximation 
such that p(AR : R) < where p is the maximum rank of the matrices describing U. 

There has recently been additional work on Robust Semidefinite Optimization, for example exploit- 
ing sparsity [115J , as well as in the area of control [751 l4"9"] . 

2.5 Robust discrete optimization 

Kouvelis and Yu [95J study robust models for some discrete optimization problems, and show that the 
robust counterparts to a number of polynomially solvable combinatorial problems are NP-hard. For 
instance, the problem of minimizing the maximum shortest path on a graph with only two scenarios for 
the cost vector can be shown to be an NP-hard problem [95] . 

Bertsimas and Sim (33] . however, present a model for cost uncertainty in which each coefficient Cj 
is allowed to vary within the interval [cj, Cj + dA, with no more than T > coefficients allowed to vary. 
They then apply this model to a number of combinatorial problems, i.e., they attempt to solve 

minimize c T x + max djXj 
{s 1 scn, \s\<t}*-' 
— — jes 

subject to x £ X, 
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where ./V = {1, . . . , n} and X is a fixed set. Under this model for uncertainty, the robust version of 
a combinatorial problem may be solved by solving no more than n + 1 instances of the underlying, 
nominal problem. This result extends to approximation algorithms for combinatorial problems. For 
network flow problems, the above model can be applied and the robust solution can be computed by 
solving a logarithmic number of nominal, network flow problems. 

Atamturk [3] shows that, under an appropriate uncertainty model for the cost vector in a mixed 
0-1 integer program, there is a tight, linear programming formulation of the robust problem with size 
polynomial in the size of a tight linear programming formulation for the nominal problem. 

3 Choosing Uncertainty Sets 

In addition to tractability, a central question in the Robust Optimization literature has been proba- 
bility guarantees on feasibility under particular distributional assumptions for the disturbance vectors. 
Specifically, what does robust feasibility imply about probability of feasibility, i.e., what is the smallest 
e we can find such that 

x € X(U) ^ F {fiix,^) > 0) < e, 

under (ideally mild) assumptions on a distribution for ml 

Such implications may be used as guidance for selection of a parameter representing the size of the 
uncertainty set. More generally, there are fundamental connections between distributional ambiguity, 
measures of risk, and uncertainty sets in robust optimization. In this section, we briefly discuss some 
of the connections in this vein. 

3.1 Probability Guarantees 

Probabilistic constraints, often called chance constraints in the literature, have a long history in stochas- 
tic optimization. Many approaches have been considered to address the computational challenges they 
pose (_ \V21\ fTTT] ). including work using sampling to approximate the chance constraints [I7J H5J [71] . 

One of the early discussions of probability guarantees in RO traces back to Ben-Tal and Nemirovski 
|15j . who propose a robust model based on ellipsoids of radius in the context of robust LP. Under this 
model, if the uncertain coefficients have bounded, symmetric support, they show that the corresponding 
robust feasible solutions must satisfy the constraint with high probability. Specifically, consider a linear 
constraint j 

< hi, where the coefficients (X% j £1X6 uncertain and given by 5y = (1 + e£jj)ajj, where 
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a>ij is a "nominal" value for the coefficient and are zero mean, independent over j, and supported 
on [—1,1]. Then a robust constraint of the form 



J \ 3 

implies the robust solution satisfies the constraint with probability at least 1 — e~^ 2 ^ 2 . This bound 
holds for any such distribution on the finite support. 

In a similar spirit, Bertsimas and Sim [34J propose an uncertainty set of the form 



(3.8) 



s||oo < i, < r 

for the coefficients a of an uncertain, linear constraint. Here, 1 : R — >• R denotes the indicator function 
of y, i.e., l(y) = if and only if y = 0, A is a vector of "nominal" values, J C {1, ... ,n} is an index 
set of uncertain coefficients, and T < |J| is an integer]^] reflecting the number of coefficients which are 
allowed to deviate from their nominal values. The dual formulation of this as a linear optimization is 
discussed in Section [2] The following then holds. 

Theorem 4. (Bertsimas and Sim Let x* satisfy the constraint 

max a T x* < b, 



where Ur is as in (3.8). If the random vector a has independent components with a,- distributed sym- 



metrically on [aj — aj, a,j + dj] if j £ J and a,j = dj otherwise, then 

f a 1 x* >b] < e'^. 



In the case of linear optimization with only partial moment information (specifically, known mean 
and covariance), Bertsimas et al. |31] prove guarantees for the general norm uncertainty model used in 
Theorem [2j For instance, when || • || is the Euclidean norm, and x* is feasible to the robust problem, 
Theorem [2] can be shown [31] to imply the guarantee 

dJx* > b] < ' 



1 + A 2 ' 

where A is the radius of the uncertainty set, and the mean and covariance are used for A and M, 
respectively. 



The authors also consider F non-integer, but we omit this straightforward extension for notational convenience. 
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For more general robust conic optimization problems, results on probability guarantees are more elu- 
sive. Bertsimas and Sim are able to prove probability guarantees for their approximate robust solutions 
in [35]. In Chen et al. [56], more general deviation measures are considered that capture distributional 
skewness, leading to improved probability guarantees. Also of interest is the work of Paschalidis and 
Kang on probability guarantees and uncertainty set selection when the entire distribution is available 



3.2 Distributional Uncertainty 

The issue of limited distributional information is central and has been the subject of considerable 
research in the decision theory literature. This work closely connects to robustness considerations and 
provides potential guidance and economic meaning to the choice of particular uncertainty sets. 

Consider a function u(x,£) where £ is a random parameter on some measure space (Q, J-). For the 
purposes of this discussion, let u be a concave, nondecreasing payoff function. In many situations, it 
may be unreasonable to expect the decision maker to have a full description of the distribution of £, but 
instead knows the distribution to be confined to some set of distributions Q. Using a well-known duality 
result that traces back to at least the robust statistics literature (e.g., Huber [89]), one can establish 
that for any set Q, there exists a convex, non-increasing, translation-invariant, positive homogeneous 
function ji on the induced space of random variables, such that 



The function in this representation falls precisely into the class of coherent risk measures popularized by 
Artzner et al. [2]. These functions provide an economic interpretation in terms of a capital requirement: 
if X is a random variable (e.g., return), n(X) represents the amount of money required to be added to 
X in order to make it "acceptable," given utility function u. The properties listed above are natural 
in a risk management setting: monotonicity states that one position that always pays off more than 
another should be deemed less risky; translation invariance means the addition of a sure amount to a 
position reduces the risk by precisely that amount; positive homogeneity means risks scale equally with 
the size of the stakes; and convexity means diversification among risky positions should be encouraged. 

The above observation implies an immediate connection between these risk management tools, dis- 
tributional ambiguity, and robust optimization. These connections have been explored in recent work 
on robust optimization. Natarajan et al. [109J investigate this connection with a focus on inferring risk 
measures from uncertainty sets. 



[116j. 




(3.9) 
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Bertsimas and Brown |23] examine the question from the opposite perspective: namely, with risk 
preferences specified by a coherent risk measure, they examine the implications for uncertainty set 
structure in robust linear optimization problems. Due to the duality above, a risk constraint of the 
form n(a!x — b) < on a linear constraint with an uncertain vector a can be equivalently expressed as 

a'x > b V a eU, 

where IA = conv ({Eq [a] : Q 6 Q}) and Q is the generating family for \i. 

For a concrete application of this, one of most famous coherent risk measures is the conditional 
value-at-risk (CVaR), defined as 



u(X) = inf \u+ -E \{-v-Xy 
v&R [ a 1 

for any a £ (0, 1]. For atomless distributions, CVaR is equivalent to the expected value of the random 
variable conditional on it being in its lower a quantile. 

Consider the case when the uncertain vector a follows a discrete distribution with support {a±, . . . , ajy} 
and corresponding probabilities {p\, ■ ■ ■ ,Pn}- The generating family for CVaR in this case is Q = {q € 
A N : qi < Pi/a}. This leads to the uncertainty set 



U = conv 



(j^& ai+ ^-^&j a i : I^{h---,N}, j G {!,..., N}\ I, X><aj^J 



This set is a polytope, and therefore the robust optimization problem in this case may be reformulated 
as a linear program. When pi = 1/N and a = j/N for some j G Z, , this has the interpretation of the 
convex hull of all j-point averages of A. 

Despite its popularity, CVaR represents only a special case of a much broader class of coherent 
risk measures that are comonotone. These risk measures satisfy the additional property that risky 
positions that "move together" in all states cannot be used to hedge one another. Extending a result 
from Dellacherie [65] . Schmeidler [124] shows that the class of such risk measures is precisely the same 
as the set of functions representable as Choquet integrals (Choquet, [59J). Choquet integrals are the 
expectation under a set function that is non-additive and are a classical approach towards dealing 
with ambiguous distributions in decision theory. Bertsimas and Brown |24j discuss how one can form 
uncertainty sets in RO with these types of risk measures on discrete event spaces. 

The use of a discrete probability space may be justified in situations when samples of the uncertainty 
are available. Delage and Ye [M] have proposed an approach to the distribution-robust problem 

minimizeajgx maxE^ [h(x, £)] , 
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where £ is a random parameter with distribution /t on some set of distributions T> supported on a 
bounded set S, h is convex in the decision variable x, and X is a convex set. They consider sets of 
distributions T> based on moment uncertainty with a particular focus on sets that have uncertainty in 
the mean and covariance of £. They then consider the problem when one has independent samples 
£i> • • • > Cm and focus largely on the set 

£>i(<S,£ , £0,71,72) = 

{p(£ G< S) = l : (E[£]-/^)'E ~>|£]-£o)<71, E [(£ - £o)'(£ " £o)] ^ 72^0} • 

The above problem can be solved in polynomial time, and, with proper choices of 71,72 and M, the 
resulting optimal value provides an upper bound on the expected cost with high probability. In the 
case of h as a piecewise linear, convex function, the resulting problem reduces to solving an SDP. This 
type of approach seems highly practical in settings (prevalent in many applications, e.g., finance) where 
samples are the only relevant information a decision maker has on the underlying distribution. 

Related to distributional uncertainty is the work in [T32]. Here, Xu, Caramanis and Mannor show 
that any robust optimization problem is equivalent to a distributionally robust problem. Using this 
equivalence to robust optimization, they show how robustness can guarantee consistency in sampled 
problems, even when the nominal sampled problem fails to be consistent. 

More general types of robust optimization models have been explored, and such approaches draw 
further connections to research in decision theory. Ben-Tal et al. [6] propose an approach called soft 



robust optimization applicable in settings of distributional ambiguity. They modify the constraint (3.9) 
and consider the more general constraint 

inf Eq [/(*,£)] > -e Ve>0, 
QeQ(e) 

where {Q(e)} e >o is a set of sets of distributions, nondecreasing and convex on e > 0. This set of 
constraints considers different sized uncertainty sets with increasingly looser feasibility requirements as 
the uncertainty size grows; as such, it provides a potentially less conservative approach to RO than 



(3.9). This approach connects to the approach of convex risk measures (Fdllmer and Schied, [H]), a 
generalization of the coherent risk measures mentioned above. Under a particular form for Q(e) based 
on relative entropy deviations, this model recovers the multiplier preferences of Hansen and Sargent 
|87j . who develop their approach from robust control ideas in order to deal with model mis-specification 
in the decision making of economic agents (see also Maccheroni et al. |1U2| for a generalization known 
as variational preferences) . 
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In short, there has been considerable work done in the domain of uncertainty set construction for 
RO. Some of this work focuses on the highly practical matter of implied probability guarantees under 
mild distributional assumptions or under a sufficiently large number of samples; other work draws 
connections to objects that have been axiomatized and developed in the decision theory literature over 
the past several decades. 

4 Robust Adaptable Optimization 

Thus far this paper has addressed optimization in the static, or one-shot case: the decision-maker 
considers a single-stage optimization problem affected by uncertainty. In this formulation, all the de- 
cisions are implemented simultaneously, and in particular, before any of the uncertainty is realized. In 
dynamic (or sequential) decision-making problems this single-shot assumption is restrictive and con- 
servative. For example, in the inventory control example we discuss below, this would correspond to 
making all ordering decisions up front, without flexibility to adapt to changing demand patterns. 

Sequential decision-making appears in a broad range of applications in many areas of engineering 
and beyond. There has been extensive work in optimal and robust control (e.g., the textbooks |67U136] , 
or the articles [721 E3 E3 E3] , and references therein) , and approximate and exact dynamic programming 
(e.g., see the textbooks \22 \ |2T | [T22j). In this section, we consider modeling approaches to incorporate 
sequential decision-making into the robust optimization framework. 

4.1 Motivation and Background 

In what follows, we refer to the static solution as the case where the x^ are all chosen at time 1 before 
any realizations of the uncertainty are revealed. The dynamic solution is the fully adaptable one, where 
Xi may have arbitrary functional dependence on past realizations of the uncertainty. 

The question as to when adaptability has value is an interesting one that has received some attention. 
The papers by Dean, Goemans and Vondrak ([62j[81]) consider the value of adaptability in the context 
of stochastic optimization problems. They show there that for the stochastic knapsack problem, the 
value of adaptability is bounded: the value of the optimal adaptive solution is no more than a constant 
factor times the value of the optimal non-adaptive solution. In [27J, Bertsimas and Goyal consider a 
two-stage mixed integer stochastic optimization problem with uncertainty in the right-hand-side. They 
show that a static robust solution approximates the fully-adaptable two-stage solution for the stochastic 
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problem to within a factor of two, as long as the uncertainty set and the underlying measure, are both 
symmetric. 

Despite the results for these cases, we would generally expect approaches that explicitly incorporate 
adaptivity to substantially outperform static approaches in multi-period problems. There are a number 
of approaches. 

Receding Horizon 

The most straightforward extension of the single-shot Robust Optimization formulation to that of 
sequential decision-making, is the so-called receding horizon approach. In this formulation, the static 
solution over all stages is computed, and the first-stage decision is implemented. At the next stage, the 
process is repeated. In the control literature this is known as open-loop feedback. While this approach 
is typically tractable, in many cases it may be far from optimal. In particular, because it is computed 
without any adaptability, the first stage decision may be overly conservative. 

Stochastic Optimization 

In Stochastic Optimization, the basic problem of interest is the so-called complete recourse problem (for 
the basic definitions and setup, see [39 1 f90l [121 j . and references therein). In this setup, the feasibility 
constraints of a single-stage Stochastic Optimization problem are relaxed and moved into the objective 
function by assuming that after the first-stage decisions are implemented and the uncertainty realized, 
the decision-maker has some recourse to ensure that the constraints are satisfied. The canonical example 
is in inventory control where in case of shortfall the decision-maker can buy inventory at a higher cost 
(possibly from a competitor) to meet demand. Then the problem becomes one of minimizing expected 
cost of the two-stage problem. If there is no complete recourse (i.e., not every first-stage decision can be 
completed to a feasible solution via second-stage actions) and furthermore the impact and cost of the 
second-stage actions are uncertain at the first stage, the problem becomes considerably more difficult. 
The feasibility constraint in particular is much more difficult to treat, since it cannot be entirely brought 
into the objective function. 

When the uncertainty is assumed to take values in a finite set of small cardinality, the two-stage 
problem is tractable, and even for larger cardinality (but still finite) uncertainty sets (called scenarios), 
large-scale linear programming techniques such as Bender's decomposition can be employed to obtain 
a tractable formulation (see, e.g., [37] ). In the case of incomplete recourse where feasibility is not 
guaranteed, robustness of the first-stage decision may require a very large number of scenarios in order 
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to capture enough of the structure of the uncertainty. In the next section, we discuss a robust adaptable 
approach called Finite Adaptability that seeks to circumvent this issue. 

Finally, even for small cardinality sets, the multi-stage complexity explodes in the number of stages 
([125J). This is a central problem of multi-stage optimization, in both the robust and the stochastic 
formulations. 

Dynamic Programming 

Sequential decision-making under uncertainty has traditionally fallen under the purview of Dynamic 
Programming, where many exact and approximate techniques have been developed - we do not review 
this work here, but rather refer the reader to the books [22], [21], and |122j . The Dynamic Programming 
framework has been extended to the robust Dynamic Programming and robust MDP setting, where 
the payoffs and the dynamics are not exactly known, in Iyengar [HI] and Nilim and El Ghaoui |113| . 
and then also in Xu and Mannor |134j . Dynamic Programming yields tractable algorithms precisely 
when the Dynamic Programming recursion does not suffer from the curse of dimensionality. As the 
papers cited above make clear, this is a fragile property of any problem, and is particularly sensitive to 
the structure of the uncertainty. Indeed, the work in [911 11131 H34[ 163] assumes a special property of 
the uncertainty set ("rect angularity") that effectively means that the decision-maker gains nothing by 
having future stage actions depend explicitly on past realizations of the uncertainty. 

This section is devoted precisely to this problem: the dependence of future actions on past realiza- 
tions of the uncertainty. 

4.2 Tractability of Robust Adaptable Optimization 

The uncertain multi-stage problem with deterministic set-based uncertainty, i.e., the robust multi-stage 
formulation, was first considered in |10j . There, the authors show that the two-stage linear problem 
with deterministic uncertainty is in general AP-hard. Consider the generic two-stage problem: 

min : c T cci 

(4.10) 

s.t. : A\{u)xx + A 2 (u)x2(u) < b, \/u£li. 

Here, X2(-) is an arbitrary function of u. We can rewrite this explicitly in terms of the feasible set for 
the first stage decision: 

min : c T x-i 

(4.11) 

s.t. : xi G {xi : Via eW, 3x 2 s.t. Ai(u)x\ + A2{u)x2 < b} . 



22 



The feasible set is convex, but nevertheless the optimization problem is in general intractable. Consider 
a simple example given in |10| : 

min : x\ 

s.t. : x\ — u T X2(u) > 

V ' ~ (4.12) 

X2(u) > Bu 
X2(u) < Bu. 

It is not hard to see that the feasible first-stage decisions are given by the set: 

{xi : x\ > u 1 Bu, Vtt G U}. 

The set is, therefore, a ray in R , but determining the left endpoint of this line requires computing a 
maximization of a (possibly indefinite) quadratic u T Bu, over the set Li. In general, this problem is 
NP-hard (see, e.g., [75]). 

4.3 Theoretical Results 

Despite the hardness result above, there has been considerable effort devoted to obtaining different 
approximations and approaches to the multi-stage optimization problem. 

4.3.1 Affine Adaptability 

In [10], the authors formulate an approximation to the general robust multi-stage optimization problem, 
which they call the Affinely Adjustable Robust Counterpart (AARC). Here, they explicitly parameterize 
the future stage decisions as affine functions of the revealed uncertainty. For the two-stage problem 



(4.10), the second stage variable, X2{u), is parameterized as: 

x 2 (u) = Qu + q. 

Now, the problem becomes: 



min : c T X\ 



s.t.: Ai(u)xi + A2(u)[Qu + q]<b, Mu^U. 

This is a single-stage RO. The decision-variables are {x\,Q,q), and they are all to be decided before 
the uncertain parameter, u £ IA, is realized. 



In the generic formulation of the two-stage problem (4.10), the functional dependence of a^O) on 



u is arbitrary. In the affine formulation, the resulting problem is a linear optimization problem with 
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uncertainty. The parameters of the problem, however, now have a quadratic dependence on the uncertain 
parameter u. Thus in general, the resulting robust linear optimization will not be tractable - consider 



again the example (4.12) 



Despite this negative result, there are some positive complexity results concerning the affine model. 
In order to present these, let us explicitly denote the dependence of the optimization parameters, A\ 
and A 2 , as: 

L 

[A u A 2 ](u) = [Af,Af ] + Y,M4\ Af]. 

1=1 

When we have A^p = 0, for all I > 1, the matrix multiplying the second stage variables is constant. This 
setting is known as the case of fixed recourse. We can now write the second stage variables explicitly 
in terms of the columns of the matrix Q. Letting denote the I th column of Q, and = q the 
constant vector, we have: 

x 2 = Qu + q 

1=1 

Letting x = ( x i,Q 1 ■ ■ ■ , Q^) denote the full decision vector, we can write the i th constraint as 

L 



< (A? ) x 1 + A¥ ) q {0) -b) l + ^u l (A^ Xl + A 2 q(% 

1=1 

L 

= Ete 

z=o 

where we have defined 

4 ± al( X ) ± (A? Xl + A«gW). a * A ^(0)^ + A (o) g(0) _ % 

Theorem 5 ([10J). Assume we have a two-stage linear optimization with fixed recourse, and with conic 
uncertainty set: 

U = {u : 3£ s.t. Viu + V 2 £ >k d} C R L , 

where K is a convex cone with dual IC* . IfU has nonempty interior, then the AARC can be reformulated 
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as the following optimization problem: 



mm : 



s.t. : -a i (x 1 ,q ( -°\...,q^)=0, % = 1, . 



, m 



V 2 X i = 



i = 1, . . . ,m 
(fX 1 + a} ) (x 1 ,q( \...,q^) > 0, i=l,...,m 
A* >k* 0, i = l,...,m. 

If the cone JC is the positive orthant, then the A ARC given above is an LP. 
The case of non- fixed recourse is more difficult because of the quadratic dependence on u. Note that 



the example in (4.12) above involves an uncertainty-affected recourse matrix. In this non-fixed recourse 



case, the robust constraints have a component that is quadratic in the uncertain parameters, U{. These 
robust constraints then become: 



Xi + 



6<0, VugW, 



which can be rewritten to emphasize the quadratic dependence on u, as 

"4V+4V 0) -&] [A^x. + A^q^ + A^q^] + [^2u kUl A^ q m 

Writing 



< 0, \/u G U. 



* r^(o) 



[Af>« 1+ 4V 0) -6]i 



/3f } (%) 



^ a [a|^) + a|V^ 



the robust constraints can now be expressed as: 



1 = 1,. ..,L 



i I, k — 1, . . . , L, 



a i {x) + 2u 1 '(3 l ( yX )+u r T i (x)u>0, V«eW. 



Theorem 6 (|10j). Let our uncertainty set be given as the intersection of ellipsoids: 



(4.13) 



U = {u : u ] (p- 2 S k )u < 1, k = 1, . . . , K}, 
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where p controls the size of the ellipsoids. Then the original AARC problem can be approximated by the 
following semidefinite optimization problem: 



Xl 






A(x) 


A(x) T 





s.t.: —\hO,i = l,...,m (4.14) 

V A(x) T 

> o, i = 1, . . . ,m 

The constant p in the definition of the uncertainty set IA can be regarded as a measure of the level 
of the uncertainty. This allows us to give a bound on the tightness of the approximation. Define the 
constant 

A 



7 

Then we have the following. 



2 In ^6 J^Rank(5jfc)Y 



Theorem 7 (|10|). Let X p denote the feasible set of the AARC with noise level p. Let Xp Pprox denote the 
feasible set of the SDP approximation to the AARC with uncertainty parameter p. Then, for 7 defined 
as above, we have the containment: 

X 1P C X^ pprox C X p . 

This tightness result has been improved; see |66j. 

There have been a number of applications building upon affine adaptability, in a wide array of areas: 

1. Integrated circuit design: In |104| . the affine adjustable approach is used to model the yield-loss 
optimization in chip design, where the first stage decisions are the pre-silicon design decisions, while 
the second-stage decisions represent post-silicon tuning, made after the manufacturing variability 
is realized and can then be measured. 

2. Comprehensive Robust Optimization: In [7J, the authors extend the robust static, as well as the 
affine adaptability framework, to soften the hard constraints of the optimization, and hence to 
reduce the conservativeness of robustness. At the same time, this controls the infeasibility of the 
solution even when the uncertainty is realized outside a nominal compact set. This has many 
applications, including portfolio management, and optimal control. 

3. Network flows and Traffic Management: In |112j . the authors consider the robust capacity expan- 
sion of a network flow problem that faces uncertainty in the demand, and also the travel time 
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along the links. They use the adjustable framework of [10], and they show that for the structure 
of uncertainty sets they consider, the resulting problem is tractable. In |lU7j . the authors consider 
a similar problem under transportation cost and demand uncertainty, extending the work in [112J. 

4. Chance constraints: In |57| . the authors apply a modified model of affine adaptability to the 
stochastic programming setting, and show how this can improve approximations of chance con- 
straints. In |70| . the authors formulate and propose an algorithm for the problem of two-stage 
convex chance constraints when the underlying distribution has some uncertainty (i.e., an am- 
biguous distribution). 

5. Numerous other applications have been considered, including portfolio management |46} 1129] . 
coordination in wireless networks [135J, robust control [84J, and model adaptive control. 

Additional work in affine adaptability has been done in |57j . where the authors consider modified linear 
decision rules in the context of only partial distibutional knowledge, and within that framework derive 
tractable approximations to the resulting robust problems. See also Ben-Tal et al. for a detailed 
discussion of affine decision rules in multistage optimization. Recently, |29| have given conditions under 
which affine policies are in fact optimal, and affine policies have been extended to higher order polynomial 
adaptability in [25ll28j . 

4.3.2 Finite Adaptability 

The framework of Finite Adaptability, introduced in Bertsimas and Caramanis |26j and Caramanis 
[52] . treats the discrete setting by modeling the second-stage variables, X2(ti), as piecewise constant 
functions of the uncertainty, with k pieces. One advantage of such an approach is that, due to the 
inherent finiteness of the framework, the resulting formulation can accommodate discrete variables. In 
addition, the level of adaptability can be adjusted by changing the number of pieces in the piecewise 
constant second stage variables. (For an example from circuit design where such second stage limited 
adaptability constraints are physically motivated by design considerations, see |103| I127] ). 

If the partition of the uncertainty set is fixed, then the resulting problem retains the structure of the 
original nominal problem, and the number of second stage variables grows by a factor of k. In general, 
computing the optimal partition into even two regions is NP-hard [26], however, if any one of the three 
quantities: (a) dimension of the uncertainty; (b) dimension of the decision-space; or (c) number of 
uncertain constraints, is small, then computing the optimal 2-piecewise constant second stage policy 
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can be done efficiently. One application where the dimension of the uncertainty is large, but can be 
approximated by a low-dimensional set, is weather uncertainty in air traffic flow management (see |26j). 

4.3.3 Network Design 

In Atamturk and Zhang the authors consider two-stage robust network flow and design, where the 
demand vector is uncertain. This work deals with computing the optimal second stage adaptability, and 
characterizing the first-stage feasible set of decisions. While this set is convex, solving the separation 
problem, and hence optimizing over it, can be NP-hard, even for the two-stage network flow problem. 

Given a directed graph G = (V, E), and a demand vector d G MX , where the edges are partitioned 
into first-stage and second-stage decisions, E = E\ U E% , we want to obtain an expression for the feasible 
first-stage decisions. We define some notation first. Given a set of nodes, S C V, let 5 + (S), 8~(S), 
denote the set of arcs into and out of the set S, respectively. Then, denote the set of flows on the graph 
satisfying the demand by 

V d = {x G Rf : x(5 + (i)) - aj(<T(i)) > Vt G V}. 

If the demand vector d is only known to lie in a given compact set hi C MY , then the set of flows 
satisfying every possible demand vector is given by the intersection V = Hdew ^d- ^ the edge set E is 
partitioned E = E\ U E% into first and second-stage flow variables, then the set of first-stage-feasible 
vectors is: 

V(E 1 ) = p| Pvoj El V d , 

d&A 

where Proj^Pd = {xe 1 '■ (xe 1 ,xe 2 ) G Vd}- Then we have: 

Theorem 8 ([5)- A vector xe 1 is an element ofV(E\) iff xe 1 {S + (S)) — xe 1 (8~(S)) > (s, for all 
subsets SCV such that 5 + (S) C E\, where we have defined = max{d(5) : d G LI}. 

The authors then show that for both the budget-restricted uncertainty model, IA = {d : Y^iev — 
ttq, d—h < d < d+h}, and the cardinality-restricted uncertainty model, U = {d : ^2 ieV \\dk~ di\\hf\ < 
r, d — h < d < d + h}, the separation problem for the set V(Ei) is NP-hard: 

Theorem 9 ([4j). For both classes of uncertainty sets given above, the separation problem for V (E\) is 
NP-hard for bipartite G(V, B). 

These results extend also to the framework of two-stage network design problems, where the ca- 
pacities of the edges are also part of the optimization. If the second stage network topology is totally 
ordered, or an arborescence, then the separation problem becomes tractable. 
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5 Applications of Robust Optimization 

In this section, we examine several applications approached by Robust Optimization techniques. 



5.1 Portfolio optimization 

One of the central problems in finance is how to allocate monetary resources across risky assets. This 
problem has received considerable attention from the Robust Optimization community and a wide array 
of models for robustness have been explored in the literature. 

5.1.1 Uncertainty models for return mean and covariance 

The classical work of Markowitz ( |105| 1106] ) served as the genesis for modern portfolio theory. The 
canonical problem is to allocate wealth across n risky assets with mean returns 6 M. n and return 
covariance matrix S £ S" + over a weight vector w £ W 1 . Two versions of the problem arise; first, the 
minimum variance problem, i.e., 



Here, r and a are investor-specified constants, and W represents the set of acceptable weight vectors 
(W typically contains the normalization constraint e T w = 1 and often has "no short-sales" constraints, 
i.e., Wi > 0, i = 1, . . . , n, among others). 

While this framework proposed by Markowitz revolutionized the financial world, particularly for the 
resulting insights in trading off risk (variance) and return, a fundamental drawback from the practi- 
tioner's perspective is that n and E are rarely known with complete precision. In turn, optimization 
algorithms tend to exacerbate this problem by finding solutions that are "extreme" allocations and, in 
turn, very sensitive to small perturbations in the parameter estimates. 

Robust models for the mean and covariance information are a natural way to alleviate this difficulty, 
and they have been explored by numerous researchers. Lobo and Boyd [98] propose box, ellipsoidal, 
and other uncertainty sets for \i and E. For example, the box uncertainty sets have the form 




(5.15) 



or, alternatively, the maximum return problem, i.e. 




(5.16) 



M 




= 1, . . . ,n 



S 



{s e §" | < E 4i < Ey, i 



l,...,n, j = l,...,n}. 
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In turn, with these uncertainty structures, they provide a polynomial-time cutting plane algorithm for 



solving robust variants of Problems (5.15) and (5.16), e.g., the robust minimum variance problem 



min < sup w : inf /x w > r, io £ W > . (5-17) 

Costa and Paiva [60| propose uncertainty structures of the form M. = conv (fj,^, . . . , jx k ), S = 



conv (Si, . . . , Efc) , and formulate robust counterparts of (5.15) and (5.16) as optimization problems 
over linear matrix inequalities. 

Tiitiincu and Koenig [130] focus on the case of box uncertainty sets for fi and S as well and show 



that Problem (5.17) is equivalent to the robust risk-adjusted return problem 

maxl inf \ll t w-Xw t 'Sw\ : w G W > , (5.18) 

where A > is an investor-specified risk factor. They are able to show that this is a saddle-point problem, 
and they use an algorithm of Halldorsson and Tiitiincu [86] to compute robust efficient frontiers for this 
portfolio problem. 

5.1.2 Distributional uncertainty models 

Less has been said by the Robust Optimization community about distributional uncertainty for the 
return vector in portfolio optimization, perhaps due to the popularity of the classical mean-variance 
framework of Markowitz. Nonetheless, some work has been done in this regard. Some interesting 
research on that front is that of El Ghaoui et al. [79] , who examine the problem of worst-case value-at- 
risk (VaR) over portfolios with risky returns belonging to a restricted class of probability distributions. 
The e-VaR for a portfolio w with risky returns r obeying a distribution P is the optimal value of the 
problem 

min {7 : P (7 < -r T tvJ < ej . (5.19) 
In turn, the authors in |79j approach the worst-case VaR problem, i.e., 

mm{V v (w) : w G W} , (5.20) 

where 



Vp(w) :-- 



minimize 7 
subject to supP (7 < — f T w) < e 
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(5.21) 



In particular, the authors first focus on the distributional family V with fixed mean fj, and covariance 



X! >- 0. From a tight Chebyshev bound due to Bertsimas and Popescu [32J, it was known that (5.20) is 
equivalent to the SOCP 



min|7 : n(e)\\Ti 1 ^ 2 w\\2 — fx 7 w < ^ 



where k(e) = w(l — e)/e; in [79], however, the authors also show equivalence of (5.20) to an SDP, and 



this allows them to extend to the case of uncertainty in the moment information. Specifically, when 



the supremum in (5.20) is taken over all distributions with mean and covariance known only to belong 



within U, i.e., (fx, S) G U, [79] shows the following: 



1. When U = conv Si), . . . , (fi h £;))> tnen (5.20) is SOCP-representable 



2. WhenU is a set of component- wise box constraints on /i and XI, then (5.20) is SDP-representable 



One interesting extension in |79j is restricting the distributional family to be sufficiently "close" to 
some reference probability distribution Pq. In particular, the authors show that the inclusion of an 
entropy constraint 

r (jp 
log— dP < d 
cflP 



in (5.20) still leads to an SOCP-representable problem, with k(c) modified to a new value «(e, d) (for the 
details, see [79]). Thus, imposing this smoothness condition on the distributional family only requires 
modification of the risk factor. 

Pinar and Tutuncii [118J study a distribution- free model for near- arbitrage opportunities, which 
they term robust profit opportunities. The idea is as follows: a portfolio w on risky assets with (known) 
mean /x and covariance 5] is an arbitrage opportunity if (1) fjJw > 0, (2) w T Ylw = 0, and (3) 
e T w < 0. The first condition implies an expected positive return, the second implies a guaranteed 
return (zero variance), and the final condition states that the portfolio can be formed with a negative 
initial investment (loan). 

In an efficient market, pure arbitrage opportunities cannot exist; instead, the authors seek robust 
profit opportunities at level 6, i.e., portfolios w such that 

H T w - eVw T T,w > 0, and e T x < 0. (5.22) 



The rationale for the system ( 5.22 ) is similar to the reasoning from Ben-Tal and Nemirovski [15J discussed 



earlier on approximations to chance constraints. Namely, under some assumptions on the distribution 
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(boundedness and independence across the assets), portfolios that satisfy (5.22) have a positive return 
with probability at least 1 — e~ e ' 2 l 2 . The authors in then attempt to solve the maximum-9 robust 
profit opportunity problem: 



sup 



1 6 : fi ' w - 6Vw T 'Sw > 0, e T -u;<o}. (5.23) 



They then show that (5.23) is equivalent to a convex quadratic optimization problem and, under mild 
assumptions, has a closed-form solution. 

Along this vein, Popescu [120] has considered the problem of maximizing expected utility in a 
distributional-robust way when only the mean and covariance of the distribution are known. Specifically, 
[120J shows that the problem 

min E R \u(x'R)] , (5.24) 
R~(jtt,S) 

where u is any utility function and n and S denote the mean and covariance, respectively, of the random 
return R, reduces to a three-point problem. |120| then shows how to optimize over this robust objective 



(5.24) using quadratic programming. 



5.1.3 Robust factor models 

A common practice in modeling market return dynamics is to use a so-called factor model of the form 

r = fi + V T f + e, (5.25) 

where r G W 1 is the vector of uncertain returns, /x G R n is an expected return vector, / G W 71 is a 
vector of factor returns driving the model (these are typically major stock indices or other fundamental 
economic indicators), V G ]^ mxn j s the factor loading matrix, and e S M. n is an uncertain vector of 
residual returns. 



Robust versions of (5.25) have been considered by a few authors. Goldfarb and Iyengar [82] consider 



a model with / G M(0, F) and e G A/"(0, D), then explicitly account for covariance uncertainty as: 

• D G S d = {D | D = diag(d), d t G [dj.dj } 

• v eS v = {V + W \ \\Wi\\ g <pi, i = l,...,m} 

• fj,£S m = {n Q + e | \e\i < 7i, i = l,...,n}, 
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where Wi = Wei and, for G >- 0, \\w\\ g = V w T Gw. The authors then consider various robust 



problems using this model, including robust versions of the Markowitz problems (5.15) and (5.16), 
robust Sharpe ratio problems, and robust value-at-risk problems, and show that all of these problems 
with the uncertainty model above may be formulated as SOCPs. The authors also show how to compute 
the uncertainty parameters G, pi, 7$, d i: di, using historical return data and multivariate regression 
based on a specific confidence level u>. Additionally, they show that a particular ellipsoidal uncertainty 
model for the factor covariance matrix F can be included in the robust problems and the resulting 
problem may still be formulated as an SOCP. 

El Ghaoui et al. [79] also consider the problem of robust factor models. Here, the authors show how 
to compute upper bounds on the robust worst-case VaR problem via SDP for joint uncertainty models 
in (fx, V) (ellipsoidal and matrix norm-bounded uncertainty models are considered). 

5.1.4 Multi-period robust models 

The robust portfolio models discussed heretofore have been for single-stage problems, i.e., the investor 
chooses a single portfolio w £ M. n and has no future decisions. Some efforts have been made on multi- 
stage problems. Ben-Tal et al. |llj formulate the following, L-stage portfolio problem: 

n+l 



maximize 



subject to x\ = r l r l x l r l — y\ + z\, i = 1, . . . , n, I = 1, . . . , L 

n n 

4+1 = r l n - + W- + \ + - fj,\)yl - ^(1 + 4)4, I = 1, . . . , L (5.26) 



i=l i=l 



Here, x\ is the dollar amount invested in asset i at time I (asset n + 1 is cash), r' _1 is the uncertain 
return of asset i from period I — 1 to period I, y\ (z\) is the amount of asset i to sell (buy) at the 
beginning of period /, and fd (uj) are the uncertain sell (buy) transaction costs of asset i at period I. 



Of course, (5.26) as stated is simply a linear programming problem and contains no reference to the 
uncertainty in the returns and the transaction costs. The authors note that one can take a multi-stage 
stochastic programming approach to the problem, but that such an approach may be quite difficult 
computationally. With tractability in mind, the authors propose an ellipsoidal uncertainty set model 
(based on the mean of a period's return minus a safety factor Q\ times the standard deviation of that 
period's return, similar to |118| ) for the uncertain parameters, and show how to solve a "rolling horizon" 
version of the problem via SOCP. 
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Pinar and Tiitiincii [118] explore a two-period model for their robust profit opportunity problem. 
In particular, they examine the problem 



sup inf sup# 

subject to e T x 1 = (r 1 ) 1 x° (self-financing constraint) (5-27) 

(H^x 1 - eJ{x i y^ 2 x l > 

e T x° < 0, 

where x % is the portfolio from time % to time i + 1, r 1 is the uncertain return vector for period 1, 
and (/x 2 ,S2) is the mean and covariance of the return for period 2. The tractability of (5.27) depends 
critically on U, but |118j derives a solution to the problem when U is ellipsoidal. 

5.1.5 Computational results for robust portfolios 

Most of the studies on robust portfolio optimization are corroborated by promising computational 
experiments. Here we provide a short summary, by no means exhaustive, of some of the relevant results 
in this vein. 

• Ben-Tal et al. provide results on a simulated market model, and show that their robust 
approach greatly outperforms a stochastic programming approach based on scenarios (the robust 
has a much lower observed frequency of losses, always a lower standard deviation of returns, 
and, in most cases, a higher mean return). Their robust approach also compares favorably to a 
"nominal" approach that uses expected values of the return vectors. 

• Goldfarb and Iyengar [82] perform detailed experiments on both simulated and real market data 
and compare their robust models to "classical" Markowitz portfolios. On the real market data, 
the robust portfolios did not always outperform the classical approach, but, for high values of the 
confidence parameter (i.e., larger uncertainty sets), the robust portfolios had superior performance. 

• El Ghaoui et al. [79] show that their robust portfolios significantly outperform nominal portfolios 
in terms of worst-case value-at-risk; their computations are performed on real market data. 

• Tiitiincii and Koenig [130] compute robust "efficient frontiers" using real- world market data. They 
find that the robust portfolios offer significant improvement in worst-case return versus nominal 
portfolios at the expense of a much smaller cost in expected return. 
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• Erdogan et al. [69J consider the problems of index tracking and active portfolio management and 
provide detailed numerical experiments on both. They find that the robust models of Goldfarb 
and Iyengar [82j can (a) track an index (SP500) with much fewer assets than classical approaches 
(which has implications from a transaction costs perspective) and (b) perform well versus a bench- 
mark (again, SP500) for active management. 

• Delage and Ye [64J consider a series of portfolio optimization experiments with market returns 
over a six-year horizon. They apply their method, which solves a distribution-robust problem 
with mean and covariance information based on samples (which they show can be formulated 
as an SDP) and show that this approach greatly outperforms an approach based on stochastic 
programming. 

• Ben-Tal et al. |6j apply a robust model based on the theory of convex risk measures to a real- world 
portfolio problem, and show that their approach can yield significant improvements in downside 
risk protection at little expense in total performance compared to classical methods. 

As the above list is by no means exhaustive, we refer the reader to the references therein for more 
work illustrating the computational efficacy of robust portfolio models. 

5.2 Statistics, learning, and estimation 

The process of using data to analyze or describe the parameters and behavior of a system is inherently 
uncertain, so it is no surprise that such problems have been approached from a Robust Optimization 
perspective. Here we describe some of the prominent, related work. 

5.2.1 Robust Optimization and Regularization 

Regularization has played an important role in many fields, including functional analysis, numerical 
computation, linear algebra, statistics, differential equations, to name but a few. Of interest are the 
properties of solutions to regularized problems. There have been a number of fundamental connections 
between regularization, and Robust Optimization. 

El Ghaoui and Lebret consider the problem of robust least-squares solutions to systems of over- 
determined linear equations [77J. Given an over-determined system Ax = b, where A S j^ mx « anc j 
b € M m , an ordinary least-squares problem is min \ \Ax — b\\. In [77], the authors build explicit models 

X 

to account for uncertainty for the data [A b]. The authors show that the solution to the ^-regularized 
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regression problem, is in fact the solution to a robust optimization problem. In particular, the solution 
to 



minimize|| Ax — b\\ + Py HasHl + 1, 



is also the solution to the robust problem 



min max II (A + AA)x - (b + A6)|l, 

x ||AA Ab\\ F <p 



where || • \\p is the Frobenius norm of a matrix, i.e., ||A||^ = y Tr(A A). 

This result demonstrates that "robustifying" a solution gives us regularity properties. This has 
appeared in other contexts as well, for example see [97|. Drawing motivation from the robust control 
literature, the authors then consider extensions to structured matrix uncertainty sets, looking at the 
structured robust least-squares (SRLS) problem under linear, and fractional linear uncertainty structure. 

In related work, Xu, Caramanis and Mannor |133j consider ^-regularized regression, commonly 
called Lasso, and show that this too is the solution to a robust optimization problem. Lasso has been 
studied extensively in statistics and signal processing (among other fields) due to its remarkable ability 
to recover sparsity. Recently this has attracted attention under the name of compressed sensing (see 
[551 EO]). In |133j . the authors show that the solution to 



minimize|| Ax — fc» 1 1 2 + A||a?||i, 

is also the solution to the robust problem 

min max II (A + AA)x — b\\, 

x WAAW^^P 

where || • Hoo^ is oo-norm of the 2-norm of the columns. Using this equivalence, they re-prove that Lasso 
is sparse using a new robust optimization-based explanation of this sparsity phenomenon, thus showing 
that sparsity is a consequence of robustness. 

In [131] . the authors consider robust Support Vector Machines (SVM) and show that like Lasso and 
Tikhonov-regularized regression, norm-regularized SVMs also have a hidden robustness property: their 
solutions are solutions to a (non-regularized) robust optimization problem. Using this connection, they 
prove statistical consistency of SVMs without relying on stability or VC-dimension arguments, as past 
proofs had done. Thus, this equivalence provides a concrete link between good learning properties of 
an algorithm and its robustness, and provides a new avenue for designing learning algorithms that are 
consistent and generalize well. For more on this, we refer to the book chapter on Robust Optimization 
and Machine Learning |54j . 
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5.2.2 Binary classification via linear discriminants 

Robust versions of binary classification problems are explored in several papers. The basic problem 
setup is as follows: one has a collection of data vectors associated with two classes, x and y, with 
elements of both classes belonging to W 1 . The realized data for the two classes have empirical means 
and covariances (fx x ,'S x ) and (fj, y ,Y} y ), respectively. Based on the observed data, we wish to find a 
linear decision rule for deciding, with high probability, to which class future observations belong. In 
other words, we wish to find a hyperplane T~L(a, b) = {z £ R n | a T z = fe}, with future classifications on 
new data z depending on the sign of a T z — b such that the misclassification probability is as low as 
possible. (We direct the interested reader to Chapter 12 of Ben-Tal et al. [p5] for more discussion on 
RO in classification problems). 

Lanckriet et al. [96] approach this problem first from the approach of distributional robustness. 
In particular, they assume the means and covariances are known exactly, but nothing else about the 
distribution is known. In particular, the Minimax Probability Machine (MPM) finds a separating 
hyperplane (a, b) to the problem 



max la : inf P a 1 x > b > a, inf P a 1 y < b > a } , (5.28) 
where the notation x ~ (n x , means the inf is taken with respect to all distributions with mean n x 



and covariance S x . The authors then show that (5.28) can be solved via SOCP. The authors then go 



on to show that in the case when the means and covariances themselves belong to an uncertainty set 
defined as follows 

X = {(z^, S,) I (fi x - /x^S-Vz - t4) < *A HE* - s SHf < p} , (5.29) 
y = {(tiy^y) I (tiy-n^V-^tiy-^Kv 2 , H^-E^H^p}, (5.30) 



that the problem reduces to an equivalent MPM of the form of (5.28). 

Another technique for linear classification is based on so-called Fisher discriminant analysis (FDA) 
|73j . For random variables belonging to class x or class y, respectively, and a separating hyperplane a, 
this approach attempts to maximize the Fisher discriminant ratio 

>2 



f(a,fj, x ,fi y ,'E x ,'Ey) := ^ T ■ J 1 ^ , (5.31) 



a T (£ x + £ y )a' 

where the means and covariances, as before, are denoted by (fj, x , S^) and (fi y , S^). The Fisher discrim- 
inant ratio can be thought of as a "signal-to-noise" ratio for the classifier, and the discriminant 

a nom := (E. + S^r 1 ^-^) 
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gives the maximum value of this ratio. Kim et al. [94] consider the robust Fisher linear discriminant 
problem 

maximize^ min f(a, fj, x , n y , S x , (5.32) 

where IA is a convex uncertainty set for the mean and covariance parameters. The main result is that 
if IA is a convex set, then the discriminant 



a* := '(Jij-K) 



is optimal to the Robust Fisher linear discriminant problem (5.32 ), where (/i*, /z*, £*, £*) is any optimal 
solution to the convex optimization problem: 



mm 



The result is general in the sense that no structural properties, other than convexity, are imposed on 
the uncertainty set IA. 

Other work using robust optimization for classification and learning, includes that of Shivaswamy, 
Bhattacharyya and Smola |126j where they consider SOCP approaches for handling missing and uncer- 
tain data, and also Caramanis and Mannor [53], where robust optimization is used to obtain a model 
for uncertainty in the label of the training data. 

5.2.3 Parameter estimation 

Calafiore and El Ghaoui [48j consider the problem of maximum likelihood estimation for linear models 
when there is uncertainty in the underlying mean and covariance parameters. Specifically, they consider 
the problem of estimating the mean x of an unknown parameter x with prior distribution M(x, P( A p )). 
In addition, we have an observations vector y ~ ftf(y,D(A.d)), independent of x, where the mean 
satisfies the linear model 

y = C{A c )x. (5.33) 

Given an a priori estimate of x, denoted by x s , and a realized observation y s , the problem at hand 
is to determine an estimate for x which maximizes the a posteriori probability of the event (x s ,y s ). 
When all of the other data in the problem are known, due to the fact that x and y are independent 
and normally distributed, the maximum likelihood estimate is given by 

S ML (A) = argmin||F(A)*- 5 (A)|| 2 , (5.34) 

X 
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where 



F(A) 



5(A) 



A 1 Ai A 



£>- 1 / 2 (A d )C(A c ) 
P- 1/2 (A n ) 



D- l l\A d )y s 
P- 1 / 2 (A p )x s 



The authors in [IB] consider the case with uncertainty in the underlying parameters. In particularly, 
they parameterize the uncertainty as a linear-fractional (LFT) model and consider the uncertainty set 



Ai 



A G A 



IAII < 1 



(5.35) 



where A is a linear subspace (e.g., M pXl3 ) and the norm is the spectral (maximum singular value) norm. 
The robust or worst-case maximum likelihood (WCML) problem, then, is 



minimize max ||.F(A)x — g(A) 
AeAi 



(5.36) 



One of the main results in [48J is that the WCML problem (5.36) may be solved via an SDP formulation. 
When A = M pXl? , (i.e., unstructured uncertainty) this SDP is exact; if the underlying subspace has more 
structure, however, the SDP finds an upper bound on the worst-case maximum likelihood. 

Eldar et al. [HE] consider the problem of estimating an unknown, deterministic parameter x based 
on an observed signal y. They assume the parameter and observations are related by a linear model 

y = Hx + w, 

where w is a zero-mean random vector with covariance C w . The minimum mean-squared error (MSE) 
problem is 

minE M|x - x\\ 2 ] . (5.37) 

X 

Obviously, since x is unknown, this problem cannot be directly solved. Instead, the authors assume 
some partial knowledge of x. Specifically, they assume that the parameter obeys 



IxIIt < L, 



(5.38) 



where ||x||E, = x T Tx for some known, positive definite matrix T G S™, and L > 0. The worst-case 



MSE problem then is 



min max E Mix — x\ 

x=Gy{\\x\\ T <L} L " 



min max \ x 1 (I - GH) 1 (I - GH)x + Tr(GC w G 

x=Gy {\\x\\ T <L} 



")}. (5.39) 



39 



Notice that this problem restricts to estimators which are linear in the observations. [68] then shows 



that (5.39) may be solved via SDP and, moreover, when T and C w have identical eigenvectors, that the 
problem admits a closed-form solution. The authors also extend this formulation to include uncertainty 
in the system matrix H. In particular, they show that the robust worst-case MSE problem 

min max E [lla; — tell 2 ! , (5.40) 

x=Gy {\\x\\ T <L, \\SH\\<p} L J 

where the matrix H + 5H is now used in the system model and the matrix norm used is the spectral 
norm, may also be solved via SDP. 

For other work on sparsity and statistics, and sparse covariance estimation, we refer the reader to 
recent work in |114j . [I], and [51] . 

5.3 Supply chain management 

Bertsimas and Thiele [36] consider a robust model for inventory control. They use a cardinality- 



constrained uncertainty set, as developed in Section 2.2 One main contribution of [36] is to show that 



the robust problem has an optimal policy which is of the (s^, S&) form, i.e., order an amount Sj- — x^ 
if Xk < Sfc an d order nothing otherwise, and the authors explicitly compute (sfc,S&). Note that this 
implies that the robust approach to single-station inventory control has policies which are structurally 
identical to the stochastic case, with the added advantage that probability distributions need not be 
assumed in the robust case. A further benefit shown by the authors is that tractability of the problem 
readily extends to problems with capacities and over networks, and the authors in [36] characterize the 
optimal policies in these cases as well. 

Ben-Tal et al. [§] propose an adaptable robust model, in particular an AARC for an inventory control 
problem in which the retailer has flexible commitments with the supplier. This model has adaptability 
explicitly integrated into it, but computed as an affine function of the realized demands. Thus, they 



use the affine adaptable framework of Section 4.3.1 This structure allows the authors in [9] to obtain 



an approach which is not only robust and adaptable, but also computationally tractable. The model is 
more general than the above discussion in that it allows the retailer to pre-specify order levels to the 
supplier (commitments), but then pays a piecewise linear penalty for the deviation of the actual orders 
from this initial specification. For the sake of brevity, we refer the reader to the paper for details. 

Bienstock and Ozbay [38 propose a robust model for computing basestock levels in inventory control. 
One of their uncertainty models, inspired by adversarial queueing theory, is a non-convex model with 
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"peaks" in demand, and they provide a finite algorithm based on Bender's decomposition and show 
promising computational results. 

5.4 Engineering 

Robust Optimization techniques have been applied to a wide variety of engineering problems. Many 
of the relevant references have already been provided in the individual sections above, in particular in 
Section [2] and subsections therein. In this section, we briefly mention some additional work in this area. 
For the sake of brevity, we omit most technical details and refer the reader to the relevant papers for 
more. 

Some of the many papers on robust engineering design problems are the following. 

1. Structural design. Ben-Tal and Nemirovski [12J propose a robust version of a truss topology design 
problem in which the resulting truss structures have stable performance across a family of loading 
scenarios. They derive an SDP approach to solving this robust design problem. 

2. Circuit design. Boyd et al. [H] and Patil et al. |117j consider the problem of minimizing delay 
in digital circuits when the underlying gate delays are not known exactly. They show how to 
approach such problems using geometric programming. See also (1Q4|, 11031 1 127] , already discussed 
above. 

3. Power control in wireless channels. Hsiung et al. |88j utilize a robust geometric programming 
approach to approximate the problem of minimizing the total power consumption subject to 
constraints on the outage probability between receivers and transmitters in wireless channels 
with lognormal fading. For more on applications to communication, particularly the application 
of geometric programming, we refer the reader to the monograph [SS], and the review articles 
[100, [TUT] . For applications to coordination schemes and power control in wireless channels, see 

4. Antenna design. Lorenz and Boyd [99] consider the problem of building an array antenna with 
minimum variance when the underlying array response is not known exactly. Using an ellipsoidal 
uncertainty model, they show that this problem is equivalent to an SOCP. Mutapcic et al. |108j 
consider a beamforming design problem in which the weights cannot be implemented exactly, 
but instead are known only to lie within a box constraint. They show that the resulting design 
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problem has the same structure as the underlying, nominal beamforming problem and may, in 
fact, be interpreted as a regularized version of this nominal problem. 

5. Control. Notions of robustness have been widely popular in control theory for several decades 
(see, e.g., Ba§ar and Bernhard [5], and Zhou et al. |136j ). Somewhat in contrast to this literature, 
Bertsimas and Brown |23] explicitly use recent RO techniques to develop a tractable approach to 
constrained linear-quadratic control problems. 

6. Simulation Based Optimization in Engineering. In stark contrast to many of the problems we 
have thus-far described, many engineering design problems do not have characteristics captured 
by an easily-evaluated and manipulated functional form. Instead, for many problems, the physical 
properties of a system can often only be described by numerical simulation. In |30j, Bertsimas, 
Nohadani and Teo present a framework for robust optimization in exactly this setting, and describe 
an application of their robust optimization method for electromagnetic scattering problems. 

6 Future directions 

The goal of this paper has been to survey the known landscape of the theory and applications of RO. 
Some of the unknown questions critical to the development of this field are the following: 

1. Tractability of adaptable RO. While in some very special cases, we have known, tractable ap- 
proaches to multi-stage RO, these are still quite limited, and it is fair to say that most adaptable 
RO problems currently remain intractable. The most pressing research directions in this vein, 
then, relate to tractability, so that a similarly successful theory can be developed as in single- 
stage static Robust Optimization. 

2. Characterizing the price of robustness. Some work (e.g., \34 \ I134j ) has explored the cost, in terms 
of optimality from the nominal solution, associated with robustness. These studies, however, 
have been largely empirical. Of interest are theoretical bounds to gain an understanding of when 
robustness is cheap or expensive. 

3. Further developing RO from a data-driven perspective. While some RO approaches build uncer- 
tainty sets directly from data, most of the models in the Robust Optimization literature are not 
directly connected to data. Recent work on this issue (|64j. |24j ) have started to lay a foundation 
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to this perspective. Further developing a data-driven theory of RO is interesting from a theo- 
retical perspective, and also compelling in a practical sense, as many real-world applications are 
data-rich. 
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